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Verification  and  Validation  of  the  Fractal  Attrition 

Equation 


Executive  Summary 

Under  the  auspices  of  The  Technical  Cooperation  Program  (TTCP)  Joint  Systems  & 
Analysis  (JSA)  Technical  Panel  3  (Joint  Concepts  and  Analysis)  research  program, 
researchers  in  New  Zealand  and  the  United  Kingdom  have  been  collaborating  on  the 
development  of  a  new  combat  attrition  equation.  This  equation  differs  from  the 
standard  approach  developed  by  Lanchester,  as  it  includes  spatial  resolution  explicitly 
in  its  formulation.  It  was  hoped  that  this  would  allow  for  a  better  representation  of  the 
impact  of  command  and  control  and  networked  forces  in  attrition  algorithms. 

The  present  work  examines  this  new  approach,  its  relationships  to  existing  systems  of 
Lanchester  equations  and  its  prospects  for  further  development. 

It  concludes  that  the  Fractal  Attrition  Equation  (FAE)  was  derived  using  the 
assumptions  and  constraints  of  the  MANA  simulation  and  that  it  does  produce  results 
that  are  consistent  with  that  simulation.  In  effect,  the  FAE  is  verified  as  a  metamodel 
for  MANA.  Furthermore,  the  fractal  spatial  distribution  for  each  side's  forces  is  shown 
to  result  in  a  fractal  temporal  distribution  for  casualties,  in  keeping  with  much  of  the 
available  historical  data.  Flowever,  this  novel  approach  to  incorporating  the  spatial  and 
temporal  resolution  into  a  single  force-on-force  equation  does  not  represent  a  radical 
departure  from  traditional  practice.  Lanchester  himself  proposed  a  manual  spatial  and 
temporal  discretisation,  leading  to  much  the  same  results. 

The  FAE  is  shown  to  be  identical  to  a  stochastic  Lanchester  equation  formulation. 
Lastly,  it  is  shown  that  the  FAE  equation  of  state  was  derived  using  mutually  exclusive 
constraints  on  the  relationship  between  the  spatial  and  temporal  resolution  applied  to 
the  battlefield.  A  corrected  equation  of  state  is  derived,  using  a  consistent 
interpretation  of  that  relationship  which  is  shown  to  be  consistent  with  those  obtained 
from  both  deterministic  and  stochastic  Lanchester  equation  approaches.  This  was 
accepted  at  the  TTCP  TSA  TP3  workshop  held  at  Fishermans  Bend  in  November  2005. 
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1.  Introduction 


A  theory  for  success  in  battle  has  long  been  sought,  arguably  since  the  time  of  Sun  Tze  half 
a  millennium  before  the  Christian  era.  Since  at  least  the  early  Twentieth  Century,  a 
quantifiable  and  quantitative  expression  of  such  a  theory,  or  model,  has  been  desired. 
While  no  longer  regarded  as  providing  an  accurate  prediction  on  the  precise  outcome  of  a 
battle,  such  models  are  used  by  analysts  to  provide  information  to  decision  makers,  who 
base  acquisition  decisions  involving  billions  of  dollars  on  them.  Thus  research  into  such 
mathematical  models  is  ongoing. 

Researchers  in  New  Zealand  and  the  United  Kingdom  have  been  collaborating  closely  on 
the  development  of  a  new  attrition  equation  for  warfare  [1].  This  equation  is  different  from 
the  standard  approach  developed  by  Lanchester  because  it  incorporates  spatial  patterns 
into  the  analysis.  It  is  hoped  that  this  in  turn  may  allow  for  the  better  representation  of  the 
value  of  command  and  control  and  networked  forces  in  attrition  algorithms.  They 
envisaged  that  this  work  could  lead  to  the  replacement  of  many  ideas  in  existing 
aggregated  combat  models.  This  work  has  become  part  of  a  research  program  under  The 
Technical  Cooperation  Program  (TTCP)  Joint  Systems  &  Analysis  (JSA)  Technical  Panel  3 
(Joint  Concepts  and  Analysis),  with  an  initial  report  published  in  July  2005  [1]. 

The  present  work  reviews  this  new  approach  and  its  relationships  to  existing  systems  of 
Lanchester  equations.  Consequently,  it  is  recommended  that  the  reader  be  familiar  with 
the  details  of  this  new  approach  [1]  to  fully  understand  the  work  presented  in  sections  5 
and  6  below.  It  is  considered  at  three  distinct  levels;  its  theoretical  derivation  and 
implications,  its  ability  to  describe  the  behaviour  of  the  Map  Aware  Non-uniform 
Automata  (MANA)  cellular  automata  simulation  and  its  ability  to  describe  the  available 
historical  data  describing  observed  battle  outcomes. 

Prior  to  describing  and  addressing  those  issues,  it  is  considered  necessary  to  briefly  review 
the  theory  behind  Lanchester' s  equations,  extensions  to  Lanchester' s  ideas  on  combat 
modelling  and  the  theory  justifying  the  application  of  the  metamodelling  approach  in 
combat  modelling. 


2.  Lanchester' s  Equations 

During  the  First  World  War  F.  W.  Lanchester  described  one  of  the  simplest,  and  most 
enduring,  mathematical  attrition  models  of  force-on-force  combat  [2].  Fie  proposed  two 
systems  of  equations,  depending  on  whether  the  fighting  was  "collective"  or  not. 
Collective  combat  between  a  Red  side  of  strength  R  and  a  Blue  side  of  strength  B  being 
described  by  the  equations: 


1 


DSTO-TR-1822 


—  =  -kBB(t),  R(0)  =  R0 
at 

=  ~kRR(t),  B(0)  =  B0 
at 

Which  result  in  the  equation  of  state: 

(r2-R2)  kB 
(Bl-B2)  kR 


The  quadratic  form  of  which  results  in  the  system  of  equations  known  as  the  Lanchester 
Square  Law. 

Individual  combat  on  the  other  hand  is  described  by  the  equations: 


—  =  -kBRB(t)R(t),  R(0)  =  R0 
dt 

~~  =  ~kRBR(t)B(t),  B(0)  =  B0 
dt 

With  its  equation  of  state: 

(R0-R)  _kBR 

(bo-B)  kRB 


(3) 


(4) 


The  linear  form  of  which  results  in  the  system  of  equations  known  as  the  Lanchester 
Linear  Law. 


Over  the  intervening  years  there  have  been  many  attempts  to  validate  the  use  of 
Lanchester' s  equations  to  describe  combat  outcomes  through  analysis  of  historical  data  [3]. 
These  efforts  have  largely  been  unsuccessful.  Part  of  the  reason  for  this  lack  of  success  is 
the  failure  of  subsequent  users  of  Lanchester's  theory  to  understand  its  inherent 
assumptions,  constraints  and  limitations.  Typically,  most  users  assume  that  the  fighting 
will  be  entirely  collective  and  that  both  sides  will  attrite  each  other  evenly,  with  strengths 
asymptotically  approaching  zero.  Yet  random  chance  is  a  factor  in  warfare  whose  effect 
can  be  modelled  by  treating  attrition  as  a  stochastic  process.  This  will  be  covered  in  greater 
detail  below.  Moreover,  the  Lanchester  Square  Law  equations  are  unstable  and  any 
random  deviation  from  its  expectation  value  will  be  amplified.  The  Square  Law  system  of 
equations  for  force  A  can  be  reduced  to  a  second  order  differential  equation  of  the  form: 


d2A{t ) 
dt 2 


cA(t)  =  0 


(5) 
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the  homogenous  solution  of  which  has  a  positive  discriminant  and  real  roots,  contrary  to 
the  claim  made  in  [4],  However,  the  conclusion  reached  in  [4],  that  the  solution  to 
Lanchester's  equation  is  unstable,  still  applies  although  the  correct  explanation  of  that 
instability  lies  in  the  discrete  nature  of  the  attrition  process  and  the  limitations  in 
describing  that  process  with  continuous  variables  [4],  As  a  result  of  this  instability,  real 
combat  outcomes  can  and  should  deviate  considerably  from  that  predicted  by  simple 
application  of  Lanchester's  equations. 

Equally,  the  assumption  of  "collective"  combat  is  unlikely  to  apply  throughout  the  entire 
battle  and  hence  actual  attrition  results  from  a  combination  of  collective  and  individual 
combats.  This  possibility  has  long  been  recognised  and  produced  many  attempts  to 
generalise  Lanchester's  system  of  equations  to  better  represent  actual  combat  results.  This 
will  also  be  covered  in  greater  detail  below,  as  will  discussion  of  the  assumptions  that 
underlay  Lanchester's  equations. 

A  number  of  these  issues  were  recognised  by  Lanchester  himself,  and  covered  in  his 
writings,  through  his  attempt  to  validate  his  system  of  equations  by  analysis  of  the  battle 
of  Trafalgar.  That  analysis  has  recently  been  revisited  [4]  and  is  worth  reviewing  here  as  it 
illustrates  a  number  of  problems  encountered  in  appropriate  use  of  Lanchester's 
equations. 

2.1  Validation 

While  real  battles  are  non-uniform  in  space  and  time,  neither  parameter  appears  in 
Lanchester's  equations.  This  results  from  an  implicit  assumption  that  all  the  forces  of  each 
side  are  able  to  engage  the  other  (or  in  other  words,  the  ability  to  engage  is  not  target 
limited)  in  the  derivation  of  the  square  law. 

However,  from  a  study  of  the  evolution  of  the  Battle  of  Trafalgar,  Lanchester  realised  that 
a  battle  is  more  properly  viewed  as  a  series  of  concurrent  and  consecutive  sub-battles 
separated  by  space  and  time.  Lanchester  proposed  that  each  of  these  sub-battles  be 
described  using  his  square  law  and  that  the  overall  casualties  of  the  battle  are  obtained  by 
summation  of  those  losses.  A  recent  review  [4]  of  this  work  has  shown  that  each  of  those 
sub-battles  can  also  be  deconstructed  further  into  a  set  of  smaller  sub-battles,  and  the 
outcome  of  the  battle  can  be  shown  to  be  strongly  dependent  on  the  ability  of  commanders 
to  fragment  the  battle  to  their  advantage. 

In  summary,  the  successful  application  of  Lanchester's  equations  to  describe  actual  battle 
results  first  requires  an  understanding  of  the  structure  of  that  battle  to  enable  a 
hierarchical  decomposition  of  the  battle  into  smaller  sub-battles  until  the  assumptions 
underpinning  Lanchester's  equations  are  met  by  those  sub-battles.  It  could  be  said  that  this 
decomposition  results  from  the  application  of  Command  and  Control  on  the  battle's 
evolution.  At  any  rate,  the  outcomes  of  those  sub-battles  can  then  be  determined  and  the 
results  propagated  into  the  starting  conditions  for  the  following  group  of  sub-battles.  This 
process  is  repeated  until  the  end  of  the  battle.  This  is  the  process  that  Lanchester 
envisaged,  but  is  not  how  most  practitioners  use  his  theory.  While  this  approach  may  do 
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better  at  reproducing  the  results  of  known  battles,  the  appropriate  deconstruction  of  a 
battle  into  its  sub-battles  for  battles  yet  to  be  fought  is  clearly  open  to  interpretation. 

2.2  Interpretation 

In  the  derivation  of  the  square  and  linear  laws,  Lanchester  interpreted  their  functional 
difference  as  resulting  from  changes  in  technology  arising  from  the  industrial  revolution. 
His  square  law  was  considered  to  describe  "modern"  combat  while  the  linear  described 
"ancient"  combat. 

This  is  not  the  interpretation  taken  by  many  subsequent  practitioners  who  subscribe  to  the 
view  that  the  linear  law  describes  attrition  resulting  from  indirect  or  area  fire,  while  the 
square  law  describes  attrition  resulting  from  direct  or  point  fire. 

However,  the  work  of  Karr  [5]  shows  that  both  the  linear  and  square  laws  are  consistent 
with  both  the  point  and  area  fire  interpretations,  which  are  thus  seen  to  represent  sufficient 
but  not  necessary  conditions  for  the  derivation  of  the  loss  rate  terms  in  equations  (1)  and 
(3).  This  point  is  so  frequently  misunderstood  by  current  practitioners  that  it  is  worth 
repeating  the  description  here  in  full  [5]. 

1.  Square  Law 

•  Point  Fire  Interpretation  Targets  are  sufficiently  numerous  or  the  ability  to 
locate  them  is  sufficiently  good  that  each  attacker  locates  targets  at  a  constant  rate. 

•  Area  Fire  interpretation  Each  attacker  engages  all  targets  in  a  certain  area  per 

unit  time  while,  targets  are  dispersed  over  a  region  maintaining  a  constant  density 
so  that  a  reduction  in  the  number  of  targets  reduces  the  area  occupied. 

2.  Linear  Law 

•  Point  Fire  Interpretation  Targets  are  sufficiently  few  or  sufficiently  difficult  to 
locate  and  attack  so  that  each  attacker  locates  targets  at  a  rate  proportional  to  the 
number  of  targets  present. 

•  Area  Fire  interpretation  Each  attacker  engages  all  targets  in  a  certain  area  per 

unit  time,  while  targets  are  dispersed  over  a  region  maintaining  a  constant  area 
occupied  so  that  a  reduction  in  the  number  of  targets  reduces  the  target  density. 

Karr  [5]  also  shows  that  the  appropriate  assumptions  that  underlay  the  derivation  of 
Lanchester's  equations  are  most  easily  understood  in  terms  of  a  stochastic  Markov  Process 
and  will  be  considered  as  part  of  that  discussion. 

2.3  Heterogeneous  Combat  Models 

The  Lanchester  equations  above  use  a  single  force  parameter  to  describe  each  side's 
strength.  This  was  a  reasonable  description  of  the  forces  of  Lanchester's  day.  It  no  longer 
holds  true  today,  in  which  each  side  contains  a  number  of  different  and  effective  force 
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elements.  A  number  of  approaches  to  deal  with  this  increased  battlefield  complexity  have 
been  developed. 

Heterogeneous  Lanchester  equations  extend  equations  (1)  and  (3)  to  systems  of  equations 
where  each  side's  single  strength  is  replaced  by  a  series  of  strengths  describing  the 
numbers  of  each  type  of  combatant.  In  this  case,  the  heterogeneous  square  law  equations 
can  be  written  as: 

N 

YJkB(i,j)Bi(t) 

(6) 

^  kR 
;=1 

In  general,  a  state  equation  no  longer  exists.  Similar  equations  can  be  developed  for  the 
Linear  Law.  This  represents  a  considerable  increase  in  complexity  of  the  description  of  the 
attrition  process  and  a  consequent  reduction  in  the  model's  utility.  As  such  we  will  not 
consider  the  heterogeneous  systems  of  equations  further. 

Most  modern  approaches  replace  this  heterogeneous  system  of  equations  with  an 
homogenous  system  of  equations  in  which  overall  force  strength  is  measured  using  some 
form  of  force  scoring.  Force  scoring  is  a  means  of  aggregating  a  side's  combat  strength 
including  both  the  number  of  weapons  of  a  given  type  and  an  assessment  of  the  relative 
contribution  that  each  weapon  type  makes  to  that  side's  combat  power.  Such  approaches 
have  been  used  with  varying  degrees  of  success  for  many  years  and  a  number  of  different 
force  scoring  methods  have  been  developed  [6]  [7], 

It  is  interesting  to  note  that  one  of  these  simple  deterministic  analyses  [8]  [9]  has  proven 
qualitatively  better  than  many  of  the  more  sophisticated  analyses  at  predicting  the 
outcomes  of  recent  conventional  armed  conflicts.  This  success  arises  from  inclusion  of  the 
effect  of  many  factors  not  considered  in  the  course  of  development  of  the  Lanchester 
Equations.  Those  effects  were  determined  from  extensive  analysis  of  historical  combat 
outcomes  and  are  frequently  derided  as  mere  curve  fitting.  Yet  the  model  is  simple  to  use 
and  at  least  as  effective  as  any  other  force  level  combat  model. 


dRt 

dt 

dBj 

dt 


3.  Attrition  Metamodels 


At  its  simplest,  a  metamodel  is  a  model  that  describes  a  model.  This  definition  is  not  as 
self -referential  as  it  might  first  appear. 

A  model  is  a  representation  of  an  actual  situation  that  may  be  used  to  better  understand 
that  situation.  Complex  phenomena  often  require  complex  models  if  the  models  behaviour 
is  to  reproduce  that  of  the  real  world.  The  modelling  of  combat  attrition  is  one  such 
situation,  as  described  above  [6]  [8],  However,  while  such  models  produce  reasonable 
agreement  with  real  world  results,  they  are  less  useful  in  understanding  the  functional 
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dependence  of  the  modelled  quantity  on  the  input  parameters.  In  such  cases  it  is  useful  to 
develop  a  (simpler)  model  of  that  model  which,  although  providing  lower  fidelity  results, 
is  better  at  explaining  the  causes  of  those  results.  Such  models  are  termed  metamodels. 

Metamodels  are  often  developed  using  an  ad-hoc  unstructured  approach  such  as 
dimensional  analysis.  In  recent  years  a  systematic  process  for  the  development  of  such 
metamodels  has  emerged  that  provides  a  degree  of  rigour  to  the  undertaking  [10]. 

While  the  interested  reader  is  referred  to  the  work  by  Barenblatt  [10]  for  a  comprehensive 
treatment  of  similarity  laws,  a  brief  summary  is  given  below.  It  must  be  noted  that  the  use 
of  dimensional  analysis  in  the  study  of  similarity  laws  is  only  strictly  correct  for  self¬ 
similar  problems.  However,  the  distribution  of  casualties  from  historical  battles  have  been 
shown  to  exhibit  self-similar  behaviour  [11],  and  hence  provides  some  justification  in  the 
application  of  this  approach  in  combat  attrition  modelling. 

The  relationship  between  an  outcome  of  the  model,  and  a  set  of  input  variables  can  be 
written  as: 

a  =  f(al,...ak,bl,..Jjm )  (7) 

where  the  variables  ai,...ak  have  independent  dimensions  (the  dimension  of  any  of  the  a's 
cannot  be  expressed  as  a  combination  of  the  dimensions  of  the  other  a's),  and  the 
dimensions  of  the  bi,...bm  can  be  expressed  in  terms  of  products  of  powers  of  the 
dimensions  of  ai,...ak.  In  general  k  >  0  and  m  >  0. 

This  allows  the  values  of  the  ai,...ak  to  be  independently  varied,  so  that 


The  dimensions  of  a  and  bi,...bm  may  then  be  represented  as  power  monomials  in  the 
dimensions  of  al , . .  ,ak ,  for  example: 


[bJ\=[a1]Pj...[ak]rj 
[a]  =  [a, 


(9) 


With  the  corresponding  transformation  of  values: 

b'rA?...A?bj 
a'=  Axp  ...A/  a 


(10) 


Introducing  the  similarity  parameters: 
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n>  =“p7 


n  = 


(ii) 


where  the  exponents  of  the  variables  are  chosen  so  that  the  parameters  n  and  II/,... II,,,  are 
dimensionless,  enables  equation  (7)  to  be  re-written  as: 


n  = 


af  ...a k 


7f(al,..Mk,nia(’\..ar^...nmalPm...arl;") 


(12) 


Barenblatt  [10]  shows  that  this  in  turn  can  be  re-written  in  terms  of  a  function  of  a  smaller 
number  of  dimensionless  variables,  leading  to  the  relationship: 


a  =  ap  ...a[  ® 


bi 


ap'  an 
Vwi  •••“* 


a,  ...a 


k  J 


(13) 


Furthermore,  self-similar  solutions  correspond  to  cases  where  the  values  of  the  variables 
bi,...bm  tend  to  zero  or  infinity.  Barenblatt  [10]  considers  three  cases: 

1.  Type  1  metamodel. 

The  function  ®  tends  to  a  non-zero  finite  limit  as  Tlj  tends  to  zero  or  infinity.  In  practice 
this  means  ®  can  be  replaced  by  its  limiting  expression,  and  hence  /  will  be  a  product  of 
power  monomials  whose  values  can  be  determined  by  dimensional  analysis. 

2.  Type  2  metamodel. 


The  function  ®  tends  to  the  power  law  asymptotic  expression: 


( 


®  =  n;® 

V 


n_ 

Ua/ 


n. 


n“ 


X 

) 


(14) 


as  Tlj  tends  to  zero  or  infinity.  The  power  law  form  of  the  limiting  expression  still  leads  to 
complete  separation  of  variables,  but  with  characteristic  exponents  which,  in  contrast  to 
the  type  1  metamodels,  cannot  all  be  determined  by  dimensional  analysis. 

3.  Type  3  metamodel. 

Power  type  asymptotic  behaviour  is  not  observed,  and  the  function  ®  has  no  finite  limit 
different  from  zero. 
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4.  Generalised  Lanchester  Equations 


As  observed  above,  Lanchester's  equations  have  generally  failed  to  correspond  with 
available  historical  data  [3].  More  limited  examination  of  the  historical  data  has  found 
related  sub-sets  of  that  data  (such  as  the  Ardennes  Campaign  and  the  Iwo  Jima  Campaign) 
that  are  consistent  with  the  different  forms  of  Lanchester's  laws  [3].  This  has  lead  many 
researchers  to  view  Lanchester's  original  equations  as  special  cases  of  a  more  general 
attrition  equation,  or  as  attrition  metamodels  as  discussed  in  the  previous  section.  Several 
different  means  for  obtaining  that  more  general  attrition  equation,  or  Generalised 
Lanchester  equation,  have  been  proposed.  They  fall  into  two  basic  categories,  depending 
on  their  fundamental  paradigm. 

4.1  Deterministic  Approaches 

There  are  two  distinct  approaches  for  deterministic  generalised  Lanchester  equations: 
arbitrary  power  law  formulations  and  non-monotonic  formulations. 

Non-monotonic  formulations  replace  the  Lanchester  equations  of  equation  (1)  with 
equations  of  the  form: 


-  =  A-kBB(t)-qRR(t),  R(0)  =  R0 
at 

~~  =  B  —  kRR(t)  —  q  BB(t),  B(0)  =  B0 
at 


(15) 


where  the  additional  terms  can  represent  reinforcements  and  decision  thresholds,  such  as 
in  the  Dewar  model  [12],  or  production  as  in  the  strategic  conflict  models  described  by 
Kimball  [13].  These  models  can  result  in  highly  non-monotonic,  even  chaotic,  behaviour 
[8]  that  does  describe  some  features  of  combat  attrition  not  seen  in  Lanchester's  laws. 
However,  they  are  in  general  no  more  successful  in  replicating  the  broad  range  of 
available  historical  data  than  the  original  Lanchester's  equations. 

The  more  common  approach  is  to  regard  Lanchester's  equations  as  metamodels  resulting 
from  the  application  of  equation  (13)  above  to  a  generalised  Lanchester  equation,  where 
the  assumptions  and  constraints  applied  to  the  generalised  equation  have  specified  the 
values  of  the  monomial  exponents.  The  general  expressions  then  have  the  form: 


=  -kcBB{t)DR{t)E,  R(0)  =  R0 
dt 

~~  =  ~kBR(t)G  B(t)H ,  B(0)  =  B0 
dt 


(16) 


where  the  exponents  ( C..H)  are  not  constrained  to  integer  values,  but  usually  D-G  =  H-E 
and  G+H  =  D+E  =  1.  This  general  expression  does  not  have  to  be  symmetric  in  B(t)  and 


8 


DSTO-TR-1822 


R(t).  It  is  also  interesting  to  note  that  the  same  equation  of  state  results  from  different 
values  of  the  exponents  if  these  two  constraints  are  met. 

Substantial  work  has  been  done  using  this  form  of  the  generalised  Lanchester  equations  by 
a  number  of  researchers  including  Hartley  [3]  and  Helmbold  [14],  to  which  the  interested 
reader  is  referred.  If  D-G  =  H-E  =  a-1  and  C-F  =  [3,  an  equation  of  state  can  be  found  for  the 
system  of  equation  (16)  such  that  [14]: 


In 

fRo-m2^ 

=  a  In 

U  J 

In [HelmboldRatio)  =  a  In [ForceRatio)+  (3 


(17) 


Hartley  [3]  has  analysed  the  results  of  over  600  battles  and  graphed  the  parameters  in  the 
equation  of  state,  as  shown  in  Figure  1. 


-3-2-10  1  2  3  4 


ln(Force  Ratio) 

VICTOR  •  •  ■  Attacker  •  •  •  Defender  - a=  1.35 

Figure  1:  Helmbold's  Relationship,  from  Hartley  [3], 
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The  line  of  best  fit  indicates  the  relationship  between  initial  and  final  strengths  that  would 
apply  if  Lanchester  like  attrition  was  the  sole  process  at  work.  The  variation  of  outcomes 
with  respect  to  that  relationship  is  due  to  variation  in  the  Lanchester  attrition  coefficients 
with  factors  such  as  cohesion,  leadership,  terrain  and  weather  that  are  not  explicitly  dealt 
with  in  his  equations.  There  have  been  a  number  of  attempts  to  determine  the  coefficient's 
dependence  on  those  properties  such  as  the  Quantified  Judgement  Model  (QJM)  [6]  and 
the  Oak  Ridge  Spreadsheet  Model  (ORBSM)  [3].  Most  attempts,  however,  have  focussed 
on  a-priori  calculation  of  the  coefficient  from  first  principles  using  the  flawed 
interpretations  of  the  form  of  Lanchester's  laws  in  representing  point  and  area  fire  already 
referred  to. 

4.2  Stochastic  Approaches 

A  continuous  time,  discrete  space,  Markov  process  can  lead  to  stochastic  attrition  models 
analogous  to  the  deterministic  Lanchester  equations  in  a  rigorous  manner.  This  has  been 
comprehensively  explored  by  a  number  of  workers  and  is  readily  accessible  [13]  [15].  The 
treatment  summarised  here  follows  that  of  Karr  [5]. 

If  X  is  a  regular  continuous  time  Markov  process  with  countable  space  E  =  N2  the  set  of 
integer  pairs  (i,j)  with  both  (i,j)  >  0,  then: 

x  =  (x,)m 

P,(x,y)=Px{x,=y}  (18) 

l:  (x,  r)-V ) 

zsE 


where  Px  is  the  probability  that  the  initial  state  (t  =  0)  is  x.  The  infinitesimal  generator 
matrix  A  is  then  defined  by: 

A(x,y)  =  l im^Lh-A.v) 

A->0  h 

p{xt+h  =  y\xt  =x}=  A(x,  y)h  +  o(h )  (19) 

p;=p,a 

A  Markov  process  (Bt,  Rt)  with  state  space  E  can  be  regarded  as  a  combat  attrition  model 
between  a  homogenous  R  side  and  a  homogenous  B  side  provided  that  the  paths  t  — >  (Bt, 
Rt)  are  non-increasing.  Here  both  Bt  and  Rt  are  random  variables  for  each  t.  Equation  (19)  is 
then  analogous  to  the  Lanchester  attrition  equations  depending  on  the  derivation  of  an 
appropriate  generator  matrix  A. 

Karr  [5]  describes  a  Markov  process  with  the  following  assumptions: 

1.  All  combatants  on  each  side  are  identical  (homogenous). 
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2.  Times  between  detections  by  a  surviving  red  combatant  are  independent  and 
identically  exponentially  distributed  with  mean  rr1,  regardless  of  the  (non-zero) 
number  of  surviving  blue  combatants. 

3.  When  a  red  combatant  detects  a  blue  combatant  an  instantaneous  attack  occurs,  in 
which  the  blue  combatant  is  destroyed  with  probability  pi  and  survives  with 
probability  1  -  pi.  Total  loss  of  contact  immediately  takes  place. 

4.  Blue  combatants  satisfy  assumptions  2  and  3  with  parameters  r2  and  pi 
respectively. 

5.  Conditioned  on  survival,  detection  and  attack  processes  of  all  combatants  are 
mutually  independent  (in  the  probabilistic  sense). 

This  can  be  described  by  the  following  generator  matrix: 

A((kj),(i~lj))  =  klj 
A((i,j),(i,j  —  1 ))  =  k2i 

,m))  =  0 . . .  V  other  pairs  of  (/,  j),  (/,  m) 

A((kj),(i,j))=~(k]j  +  k2i) 

Where  the  assumptions  above  define  the  loss  rate: 

ki=riPi...i  =  \,  2  (21) 

This  is  recognisable  as  analogous  to  Lanchester's  Square  Law.  A  similar  treatment,  albeit 
with  slightly  different  assumptions,  results  in  Lanchester's  Linear  Law  [5]. 

Karr  goes  on  to  show  that  the  expectation  values  of  the  forces'  strengths  do  not  in  general 
satisfy  the  deterministic  Lanchester  Equations  in  the  square  law  case.  However,  given  the 
additional  assumption  that  the  loss  rates  are  small,  an  equation  of  state  for  the  continuous 
time  Markov  process  can  eventually  be  found  from: 

jE[kXBPB)-k,(RPR^=  0  (22) 


5.  Fractal  Attrition  Equations 

Lauren's  combat  attrition  metamodel  [1]  was  developed  from  a  description  of  the  attrition 
process  in  the  MANA  cellular  automata  combat  model.  The  key  process  underlying  the 
model  is  the  means  by  which  agents  from  either  side  discover  opposition  agents.  This  is  a 
probabilistic  process  and  hence  results  in  a  discrete  time  stochastic  model. 


During  a  time  step  of  duration  At,  an  agent  travelling  at  a  velocity  v  will  travel  a  distance  d 
=  vAt.  The  agent  is  then  able  to  search  an  area,  for  simplicity  assumed  to  be  a  square  box 
with  dimension  d,  for  opposition  agents  during  that  time  interval.  If  the  opposition  force  is 
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treated  as  consisting  of  clusters  of  agents,  and  the  clusters  have  a  fractal  spatial 
distribution  described  by  a  fractal  dimension  D,  then  the  probability  that  at  least  one 
opposition  agent  will  be  encountered  in  that  box  is  proportional  to  d~D.  If  the  average  size 
of  a  cluster  is  n d,  then  the  average  number  of  opposition  agents  encountered  is 
proportional  to  d~Dnd.  If  R(t)  is  the  number  of  opposition  agents  and  (p(x,t)  is  the  number  of 
clusters,  then  the  expectation  value  for  the  number  of  agents  killed  by  opposition  agents 
during  that  time  step  is: 


"A B'' 

v  A t  j 


oc 


k{Atf  D  R(t)cp(x,t) 


(23) 


Invoking  Barenblatt's  metamodelling  approach  [10]  above,  this  can  be  re-written  as: 


A B 
At 


f(k,{Atf-D  ,R(t)q>{x,t)) 


(24) 


Relaxing  the  assumptions  that  were  used  to  obtain  this  general  attrition  equation,  applying 
Barenblatt's  asymptotic  expansion  process  and  constraining  the  values  of  the  exponents 
with  dimensional  analysis,  leads  to  the  result: 


A  A 


" 0<l<T 


At 


oc  kaAta-lR(t) 


(25) 


where  a  =  D/2.  The  case  a  =  1  can  be  shown  to  correspond  with  the  assumptions  of 
Lanchester's  Square  Law.  Equation  (25)  then  reduces  to  the  expected  Lanchester  result. 

This  model  is  usually  specified  in  terms  of  the  Fractal  Attrition  Equation  (FAE): 


-‘0<t<T 


/  AB^ 

f  T  -a-T>  ^ 

-  f 

—  ^0  <t<T 

-ck2  At  2  R(t) 

V  At ) 

K  ) 

(26) 


5.1  Interpretation 

If  the  combat  system  is  assumed  to  be  self-similar,  then  multifractal  statistics  apply.  This 
means  that  the  form  of  expression  above  applies  in  an  analogous  way  to  all  of  the 
moments  of  the  random  variable.  Thus  we  have: 


E 


0  <t<T 


B(t  +  At)-B(t) 
At 


)  oc  A tsiD’p) 


(27) 


In  particular,  consideration  of  the  Fourier  transform  of  the  second  moment  of  the  random 
variable  yields  a  power  spectrum  for  the  frequency  distribution  of  a  side's  casualties: 
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|3(sX/)2h  l/r'"'1  (28) 

In  general,  the  precise  form  of  the  power  law  exponent  has  to  be  determined 
experimentally.  Examination  of  available  historical  data  is  consistent  with  this  relationship 
[16]- 

5.2  Fractal  Dimension 

The  introduction  of  a  fractal  dimension  to  describe  the  spatial  distribution  of  a  side's  units 
in  equation  (23)  was  empirical.  That  is,  there  was  no  a-priori  quantification  of  how  that 
distribution  influenced  the  attrition  process,  just  the  qualitative  observation  that  the 
available  historical  data  was  consistent  with  that  assumption.  It  remains,  therefore,  to 
describe  how  a  fractal  dimension  can  be  specified  that  describes  the  differences  between 
different  spatial  distributions  of  units  in  combat  situations. 

A  number  of  different  techniques  for  the  measurement  of  a  fractal  dimension  have  been 
developed  [17],  It  is  by  no  means  clear  which  is  the  more  appropriate  or  accurate. 
However,  given  the  close  relationship  between  the  FAE  metamodel  and  the  MANA 
cellular  automata  simulation,  a  box  counting  technique  that  is  appropriate  to  MANA  was 
chosen  [1]: 


D  =  lim 


d 


In  A 

ln(l  Td) 


(29) 


where  N  is  the  number  of  boxes  necessary  to  cover  all  the  automata. 

5.3  Validation 


MANA  [11]  is  different  from  traditional  cellular  automata  models  in  that  it  includes 
limited  global  interactions  between  automata  as  well  as  local  interactions.  The  automata's 
behaviours  are  governed  by  sets  of  parameters  that  determine  whether  they  move  towards 
friendly  or  enemy  units,  towards  an  objective  and  their  interactions  with  other  automata. 
Additional  parameters  allow  event  driven  conditional  modifiers  to  these  behaviours,  such 
as  requiring  certain  conditions  to  be  met  before  an  automaton  can  move  towards  an 
objective  or  foe.  A  final  set  of  parameters  describes  the  basic  capabilities  of  the  automaton, 
such  as  weapons  range,  detection  range  and  movement  rate. 

Global  interactions  are  facilitated  by  providing  each  side  a  map  on  which  the  locations  of 
enemy  automata  can  be  stored.  The  automata  react  to  both  units  that  are  detected  directly 
and  those  they  are  aware  of  on  the  map.  The  map  thus  simulates  a  rudimentary 
communications  network. 


MANA  automata  use  a  penalty  function  to  rank  possible  moves  based  on  their  rules  set.  If 
several  moves  have  a  similarly  low  penalty  value,  one  of  them  is  chosen  at  random.  A 
parameter  sets  the  size  of  the  margin  for  discriminating  between  similarly  valued  moves. 
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Interactions  between  automata  allow  complex  combat  behaviours  to  be  simulated, 
including: 

1.  The  automata  of  each  side  are  attracted  towards  a  goal  behind  its  opponent. 

2.  Initially,  neither  side  is  in  contact  with  the  other. 

3.  No  automaton  may  advance  without  the  presence  of  at  least  some  friendly 
automatons  within  detection  range. 

4.  Once  enemy  automata  are  detected,  a  minimum  number  of  friendly  automata  is 
required  within  detection  range  to  continue  advancing. 

5.  An  automaton  that  is  unable  to  advance  will  retreat. 

6.  Once  no  enemy  are  within  detection  range,  retreating  automata  may  regroup  and 
advance  again. 

Hence  MANA  is  a  discrete  time  stepped  stochastic  simulation  that  allows  concentration 
and  reorganisation  of  units  in  both  space  and  time.  Attrition  results  from  consecutive 
stochastic  detection  and  engagement  processes.  The  continuous  reorganisation  of  the 
forces  dispositions  has  similarities  with  withdrawal  and  reinforcement  of  combat  forces. 
The  processes  that  underlie  this  reorganisation  result  from  a  combination  of  local 
knowledge  and  Command  and  Control  (C2)  effects  from  use  of  the  map  mechanism.  The 
relationship  between  the  MANA  simulation's  design  and  the  assumptions  made  in 
developing  the  FAE  should  be  immediately  apparent. 

Comparing  the  results  of  MANA  simulations  of  a  number  of  scenarios  with  the 
implications  of  the  FAE  metamodel  has  been  undertaken  [1],  Comparing  a  metamodel 
with  a  simulation  is  straightforward  as,  in  principle,  all  the  parameters  in  the  simulation 
can  be  systematically  varied  (such  as  k).  This  enables  the  full  parameter  space  to  be 
explored  for  consistency  with  the  metamodel.  These  comparisons  have  suggested  that,  as  a 
general  rule,  Lanchester's  Square  Law  was  more  accurate  than  the  FAE  metamodel  when 
the  automata  maintained  a  single  clump-like  formation,  or  when  there  was  a  disparity  in  k 
between  the  two  sides.  In  all  other  cases  the  FAE  metamodel  was  more  accurate. 

Such  comparisons  with  historical  data  are  much  more  difficult.  For  example,  k  is  then  an 
artefact  of  the  experimental  data  chosen  and  cannot  be  independently  varied. 

Beginning  with  the  FAE,  the  loss  exchange  ratio  (LER)  for  both  sides  can  be  obtained: 

(B(t  +  At)  -  B(t))  =  (kRAt)D«'2R(t) 

(R(t  +  At)-R(t))  (kBAt)Ds/2B(t) 

Given  the  assumption  that  the  entities  have  a  fractal  distribution  with  cell  size  d  being 
chosen  so  that  each  cell  contains  at  most  one  entity,  then: 

Rcc-^^D  =  c-\ogdR  (31) 

d 
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is  in  general  true  and  can  thus  be  evaluated  at  time  t  =  0.  Further,  if  the  constants  relating 
to  the  red  and  blue  forces  are  similar,  the  equation  of  state: 


AS 

oc 

w 

M 

UJ 

Q 


(32) 


can  be  obtained  and  used  in  comparison  with  the  available  historical  data.  The  exponent  Q 
is  as  defined  in  [1], 

6.  The  Fractal  Attrition  Equation;  a  Better 
Representation? 

The  FAE  has  been  derived  using  the  assumptions  and  constraints  that  describe  the 
operation  of  the  MANA  simulation.  It  has  been  shown  to  produce  results  that  are 
consistent  with  those  of  the  MANA  simulation  [1],  To  all  effective  purposes,  the  FAE  has 
been  verified  as  providing  a  suitable  metamodel  for  the  MANA  simulation. 

The  assumption  of  a  fractal  spatial  distribution  for  each  side's  forces  was  shown  to  result 
in  a  fractal  temporal  distribution  for  casualties,  in  keeping  with  much  of  the  available 
historical  data.  Furthermore,  Lanchester's  Square  Law  is  shown  to  be  the  limiting  case  for 
the  FAE  on  application  of  Lanchester's  additional  assumptions. 

Reference  [1]  suggests  that  the  FAE  metamodel  could  prove  a  more  accurate 
representation  of  combat  than  the  original  Lanchester  equations,  because  it  includes 
factors  that  are  relevant  to  detection/ manoeuvre  as  well  as  to  firepower.  This  issue  is 
explored  further  in  this  section. 

6.1  Dependent  vs.  Independent  Dimensions 

Barenblatt's  metamodelling  technique  emphasises  the  difference  between  dependent  and 
independent  variables.  In  this  context,  an  independent  variable  is  a  variable  whose 
dimension  cannot  be  expressed  as  a  combination  of  the  dimensions  of  the  other 
independent  variables.  The  asymptotic  limits  of  independent  and  dependent  variables  are 
treated  quite  differently,  as  seen  in  equation  (13). 

The  FAE  model  is  a  two  step  process  involving  detection  followed  immediately  by 
combat.  The  detection  process  is  stochastic  and  involves  the  distance  d  =  vAt.  However,  in 
the  derivation  this  is  replaced  immediately  with  the  time  interval  parameter,  resulting  in 
the  generalised  attrition  equation  (24).  This  is  where  the  difficulty  occurs.  Equation  (24) 
has  two  parameters  ( k ,  At)  that  have  time  as  their  dimension.  Only  one  of  them  should  be 
treated  as  a  power  law  monomial  in  the  asymptotic  approximation,  in  accordance  with 
Barenblatt's  approach  [10]  and  equation  (13).  Equation  (24)  was  derived  after  this 
relationship  between  space  and  time  (d  =  vAt)  was  defined.  What  should  have  occurred 
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was  that  this  constraint  be  applied  after  the  asymptotic  approximation  was  made.  Hence, 
equation  (23)  should  be  replaced  with: 


V^J  =  f(k>d~D>R(tMx,t)) 


(33) 


which  has  suitable  independent  variables.  Application  of  the  intermediate  asymptotic 
technique  yields  in  place  of  equation  (25): 


J0<t<T 


A B 
At 


oc  kadfiR(t ) 


(34) 


The  relationship  between  space  and  time  intervals  can  now  be  applied,  with  the  velocity 
term  being  absorbed  into  the  proportionality  constant  and  dimensional  analysis  used  to 
constrain  the  values  of  a  and  /?.  This  results  in  the  FAE  of  equation  (26)  as  before.  Hence 
this  minor  inconsistency  in  derivation  does  not  affect  the  FAE  model's  results. 

6.2  To  Be  or  not  To  Be  Lanchester's  Equations 

Reference  [1]  argues  that  the  FAE  is  sufficiently  different  to  Lanchester's  Square  Law,  on  a 
number  of  grounds,  that  it  warrants  consideration  as  a  replacement  for  Lanchester's 
Equation.  Among  these  differences,  the  FAE  incorporates  a  measurement  of  the  effect  of 
the  spatial  distribution  of  forces  and  the  consequent  Loss  Exchange  Ratio  (LER)  exhibits 
different  behaviours  in  time.  However,  reference  [1]  also  shows  that  Lanchester's  Square 
Law  is  the  limiting  case  for  the  FAE,  given  a  fractal  dimension  of  2  corresponding  to  a 
uniform  force  distribution.,  and  suggests  that  the  fractal  dimension  of  the  spatial 
distribution  of  forces  might  be  a  measure  of  the  state  of  its  Command,  Control, 
Communications,  Computers,  Intelligence,  Surveillance  and  Reconnaissance  (C4ISR) 
effectiveness. 


As  previously  mentioned,  many  subsequent  users  of  Lanchester's  theory  have  failed  to 
understand  how  it  was  intended  to  be  used.  In  particular,  Lanchester  realised  that  a  battle 
is  more  properly  viewed  as  a  series  of  concurrent  and  consecutive  sub-battles  separated  by 
space  and  time.  Lanchester  proposed  that  each  of  these  sub-battles  be  described  using  his 
square  law  and  that  the  overall  casualties  of  the  battle  are  obtained  by  summation  of  those 
losses.  Hence  Lanchester  intended  that  the  spatial  and  temporal  distribution  of  forces  be 
taken  into  account  in  the  application  of  his  theory.  He  just  felt  that  it  was  appropriate  to  do 
this  manually  with  a  discrete  number  of  sub-battles  corresponding  to  identifiable  phases 
in  an  overall  engagement.  The  success  of  the  FAE  is  to  describe  a  process,  for  a  wide  range 
of  situations,  in  which  this  ensemble  averaging  of  sub-battles  can  be  analytically 
incorporated  into  the  attrition  equations  themselves.  This  is  a  significant  addition  to 
Lanchester's  theory,  but  does  not  by  itself  invalidate  that  theory.  Indeed,  it  seems  just  as 
likely  that  the  FAE  may  be  no  more  than  another  of  the  class  of  generalised  Lanchester 
Equations  corresponding  to  a  different  set  of  assumptions. 
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The  similarities  do  not  end  here.  In  the  limit  as  the  time  interval  tends  to  zero,  the  FAE  of 
equation  (26)  can  be  written  more  generally  as: 


-  =  -kBG(t)B(t),  R(0)  =  R0 
dt 

~~  =  ~kRG(t)R(t),  B(0)  =  B0 
dt 


(35) 


Where  G(t)  describes  the  time  dependence  of  the  attrition  rate  and  can  be  considered  a 
measure  of  how  the  intensity  of  combat  changes  over  time.  The  equations  have,  for 
simplicity,  been  assumed  symmetric  for  blue  and  red  forces.  In  the  FAE  this  behaviour  is 
the  reason  why  the  LER  can  remain  a  constant,  in  spite  of  changing  force  ratios,  and  in 
agreement  with  the  MANA  simulation  [1]. 

However,  Koopman  [18]  has  shown  that  the  time  scale  is  a  rather  arbitrary  parameter  in 
Lanchester's  theory  and  can  be  varied  in  accordance  with  Barenblatt's  gauge 
transformation  [10]  technique,  without  changing  the  underlaying  functional  relationships 
[19].  Therefore,  introducing  a  new  time  scale  rsuch  that 


z  =  c 


i 

J  G{t)dt 


,  dx  dx  dz 

and  —  = - 

dt  dz  dt 


which  on  substitution  into  equation  (35)  yields: 


^-  =  -^LB(z),  R(0)  =  R0 
dz  c 

=  B(0)  =  B0 

dz  c 


(36) 


(37) 


This  is  of  course  Lanchester's  Square  Law  relationship.  The  transformation  of  equation 
(36)  can  be  interpreted  as  a  "linearisation"  of  the  combat  in  time.  The  combat  then  has  a 
constant  intensity  when  measured  against  time  r,  but  varies  in  intensity  against  time  t.  The 
explicit  time  independence  of  the  attrition  process  is  one  of  Lanchester's  assumptions,  so  it 
is  useful  to  note  that  the  FAE  reduce  to  this  form  when  the  time  dependence  of  the  combat 
is  incorporated  into  the  definition  of  the  time  variable  itself  instead  of  being  left  explicit. 

However,  it  does  mean  that  the  FAE  are  a  form  of  Lanchester's  Equations. 

6.3  Stochastic  Lanchester  Equation 

The  treatment  of  Lanchester  attrition  as  a  continuous  time  Markov  process  is  described  in 
equations  (18)  to  (21)  above.  The  attrition  loss  rate  is  written  as  a  product  of  two  terms;  the 
probability  of  detection  of  an  enemy  and  the  probability  of  destruction  of  a  detected 
enemy.  The  particular  form  of  the  probability  for  detection  depends  on  the  model  behind 
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the  detection  process.  Karr  [5]  illustrates  the  process  for  an  exponentially  distributed 
probability  of  detection. 

The  FAE  was  developed  using  a  different  set  of  assumptions,  in  which  the  battlefield  of 
area  A  is  divided  into  a  number  of  boxes  (with  dimension  d)  that  are  systematically 
searched  in  turn.  This  model  is  described  by  a  geometric  probability  distribution  with  a 
constant  probability  of  a  friendly  unit  finding  an  enemy  unit  in  any  given  box  of  r.  The 
resulting  mean  time  between  detections  follows  from  the  geometric  probability 
distribution  as  1/r. 

A  Markov  attrition  process,  with  a  geometric  probability  distribution  describing  the 
detection  process,  results  in  the  same  set  of  equations  (18)  to  (21)  except  that  the  attrition 
parameters  (equation  (21))  are  obtained  using  the  geometric  probability  distribution 
values.  It  just  remains  to  determine  the  probability  for  finding  an  enemy  unit  in  any  given 
box  in  this  model. 


If  there  are  N  boxes  and  the  probability  of  occupation  for  a  box  is  p,  then  for  a  geometric 
probability  distribution,  the  mean  number  of  occupied  boxes  is  Np.  If  the  force  has  a 
fractal  spatial  distribution  described  by  a  fractal  dimension  D,  then  the  number  of 
occupied  boxes  is  given  by  d'D.  This  results  in  the  probability  that  an  enemy  unit  is 
detected  in  any  given  box  of: 


1  _  d2~D 
Nd°  ~  A 


(38) 


which,  on  remembering  that  this  model  assumes  d  =  vAt,  and  substituting  into  equation 
(21)  yields  an  expression  consistent  with  equation  (23),  the  FAE  precursor  expression. 
Hence  the  FAE  can  be  derived  as  a  stochastic  Lanchester  equation,  given  the  assumptions 
that  the  detection  process  has  a  geometric  probability  distribution  and  a  fractal  spatial 
distribution.  The  FAE,  equation  (25),  might  just  as  well  be  called  the  Fractal  Lanchester 
Equation. 

6.4  Fractal  Equation  of  State 

It  has  been  demonstrated  above,  that  for  each  of  the  attrition  expressions  corresponding  to 
the  different  forms  of  Lanchester' s  Square  Law  examined  above,  there  is  an  analogous 
equation  of  state. 

Table  1:  Model  Correspondences 


Model 

Attrition  Equation  Reference 

Equation  of  State  Reference 

Lanchester  Square  Law 

(1) 

(2) 

Generalised  Lanchester 

(16) 

(17) 

Stochastic  Lanchester 

(18) 

(22) 

FAE 

(26) 

(32) 

It  should  also  be  noted  that  when  the  assumptions  behind  the  alternate  forms  of 
Lanchester's  Square  Law  were  relaxed  to  agree  with  those  of  Lanchester's  Square  Law,  the 
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limiting  form  of  the  corresponding  attrition  expression  also  agreed  with  Lanchester's 
Square  Law. 

This  agreement  with  Lanchester's  Square  Law  as  the  limiting  case  should  also  occur  for 
the  corresponding  equations  of  state,  after  all  the  attrition  expressions  from  which  they  are 
derived  agree  in  this  limit.  Agreement  was  observed  for  the  Generalised  Lanchester  model 
and  the  Stochastic  Lanchester  model,  once  the  additional  assumption  of  large  forces  that 
was  made  in  its  derivation  [5]  is  understood. 

This  agreement  in  the  limit  is  not  observed  for  the  original  FAE  Equation  of  State  of 
reference  [1],  equation  (32).  The  Lanchester  Square  Law  Equations  of  State  are  expressions 
that  involve  the  square  of  the  forces'  strengths,  hence  the  name.  However,  the  FAE 
Equation  of  State,  resulting  from  its  Square  Law  attrition  expression,  is  linear  in  the  forces' 
strengths. 

Examining  the  derivation  of  equation  (32)  in  [1],  it  was  seen  that  the  spatial  resolution  of 
the  battlefield  d  was  chosen  as  if  it  were  independent  of  the  time  interval  At,  in  the 
derivation  of  an  approximation  expression  for  the  fractal  dimension  D.  This  is  in  contrast 
to  the  derivation  of  the  attrition  equation  (26)  where  the  relationship  between  the  time 
interval  and  spatial  resolution  was  fundamental  to  the  form  of  the  FAE.  The  consequence 
of  this  inconsistency  was  that  the  time  interval  was  varied  independently  of  the  spatial 
resolution,  and  hence  was  used  to  relate  the  initial  and  final  strength  values  in  the 
equation  of  state  resulting  in  its  linear  form. 

The  resolution  of  this  inconsistency  comes  from  recognition  that  specifying  the  spatial 
resolution,  and  hence  the  fractal  dimension,  also  determine  the  time  interval.  If  we  take  d 
to  be  the  length  of  the  cell  of  the  smallest  area  which  at  most  contains  a  single  Red  entity, 
and  if  we  assume  that  d  and  k  are  about  the  same  for  the  Red  and  Blue  forces,  substituting 
equation  31  into  the  right  hand  side  of  equation  (30)  produces: 

l^l(k  At)logAR°/B°)/2  (39) 

B(ty 


which  can  be  rearranged  to  give: 


m 

B(t ) 


f  B  ^ 

°o 

V^oy 


Q 


(40) 


where  Q  has  the  form: 

£>  =  0.5  \ogd{k  At) 


(41) 


Replacing  the  right  hand  side  of  equation  (30)  with  equation  (40),  and  cross-multiplying 
yields: 
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B(t)  AB  = 


(  n 
R0 

yRo  j 


R(t)AR 


(42) 


Dividing  the  interval  between  the  beginning  and  ending  of  the  engagement  into  a  series  of 
time  steps  allows  equation  (40)  to  be  applied  to  each  of  these  time  steps.  A  connection 
between  the  initial  and  final  force  strengths  can  then  be  obtained  by  summing  all  those 
expressions: 


Y,bab 

int  ervals 


f  n  \Q 
ZSL 

v^oy 


ZR&R 

int  ervals 


(43) 


In  keeping  with  the  development  of  the  FAE  from  a  metamodel  for  MANA,  the  force 
strengths  represent  the  number  of  live  entities.  When  attrition  occurs,  that  number  is 
reduced.  This  suggests  that  a  stochastic  model  is  more  appropriate.  If  we  now  apply  this 
additional  assumption,  that  attrition  is  modelled  as  a  Markov  attrition  process  of  the  type 
examined  by  Karr  [5],  the  probability  of  a  change  in  the  state  (force  strength)  can  be 
obtained  from  equation  (19).  The  form  of  the  infinitesimal  generator  matrix  only  allows 
transitions  to  states  whose  strengths  differ  from  the  initial  state  by  1.  The  assumption  of  a 
small  At  can  be  interpreted  as  equivalent  to  the  probability  of  more  than  one  transition  in 
any  given  time  interval  being  small,  and  which  in  the  first  instance  can  be  neglected. 
Remembering  the  connection  between  the  time  interval  and  spatial  resolution,  the  time 
interval  At  can  be  chosen  small  enough  such  that  its  corresponding  spatial  resolution  also 
satisfies  the  requirements  of  equation  (31).  Then  in  any  interval  in  the  sum  terms  of 
equation  43  in  which  there  is  no  state  transition,  either  A B  or  AR  is  0  and  equation  (43) 
reduces  to: 


I* 


B=B 


o 


f  n  \Q 
ZAL 

KRo  J 


I* 


R=R 


0 


Evaluating  the  sums,  results  in: 


f 

In 

V 


[B2{t)  +  B{t))-{Bl+BS 
f?2(0  +  i?(o)-(*o+*o)J 


f 

=  Ql  n 

v 


Bn 


R, 


o  J 


(44) 


(45) 


This  has  obvious  similarities  to  the  equations  of  state  from  both  the  Stochastic  and 
Generalised  Lanchester  models,  and  has  Lanchester's  Square  law  as  the  limiting  case  as 
required.  It  is  also  in  agreement  with  the  comprehensive  historical  data  of  Figure  1. 
Furthermore,  by  relaxing  the  assumption  that  the  value  of  k  is  about  the  same  for  the  Red 
and  Blue  forces,  and  setting  kR  =  rkB  in  equation  (30),  the  equation  of  state  can  be  re¬ 
evaluated  as: 
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In 


( B2(t)  +  B(t))-(B20+B0 ) 
(R2(t)  +  R(t))-(R>+R0) 


=  Q  In 


^  D  A 
ZA 

\R0  j 


+  0.5 Dr  In  r 


(46) 


The  similarity  between  this  equation  and  the  relationship  developed  by  Helmbold 
(equation  (17)  above)  to  describe  his  analysis  of  real-world  air  and  land  combat  results, 
and  used  by  Hartley  in  his  work  on  combat  modelling  (in  Figure  1),  is  obvious.  In  the  cases 
considered  by  Helmbold  and  Hartley,  the  forces  strengths  are  sufficiently  large  such  that 
the  difference  between  B 2  and  B(B+1)  would  not  be  observed. 

This  last  expression  is  proposed  as  the  correct  equation  of  state  for  the  FAE  model,  as  a 
replacement  for  equation  (32).  This  was  accepted  at  the  TTCP  JSA  TP3  workshop  held  at 
Fishermans  Bend  in  November  2005. 


7.  Conclusions 


The  FAE  has  been  derived  using  the  assumptions  and  constraints  that  describe  the 
operation  of  the  MANA  simulation  [1].  Furthermore,  it  has  been  shown  to  produce  results 
that  are  consistent  with  those  of  the  MANA  simulation.  For  all  effective  purposes,  the  FAE 
has  been  verified  as  providing  a  suitable  metamodel  for  the  MANA  simulation. 

The  assumption  of  a  fractal  spatial  distribution  for  each  side's  forces  was  shown  to  result 
in  a  fractal  temporal  distribution  for  casualties,  in  keeping  with  much  of  the  available 
historical  data.  Furthermore,  Lanchester's  Square  Law  was  shown  to  be  the  limiting  case 
for  the  FAE  on  application  of  Lanchester's  additional  assumptions. 

The  FAE  uses  a  novel  approach  to  incorporating  the  spatial  and  temporal  variations  in 
combat  intensity  into  a  single  force-on-force  equation.  However,  this  does  not  represent  a 
radical  departure  from  traditional  practice.  Lanchester  himself  proposed  a  paradigm  that 
involved  manually  dividing  the  battlefield  into  sub-battles  and  aggregating  their  results. 
Indeed,  it  is  easy  to  show  by  defining  a  new  time  variable  that  the  FAE  reduces  to 
Lanchester's  form. 

Furthermore,  the  FAE  can  be  derived  as  a  stochastic  Lanchester  equation,  given  the 
assumptions  that  the  detection  process  has  a  geometric  probability  distribution  and  a 
fractal  spatial  distribution.  The  FAE  might  just  as  well  be  called  the  Fractal  Lanchester 
Equation. 

Lastly,  it  is  shown  that  the  original  FAE  equation  of  state  was  derived  using  mutually 
exclusive  constraints  on  the  relationship  between  the  spatial  and  temporal  resolution 
applied  to  the  battlefield.  A  new  equation  of  state,  using  a  consistent  interpretation  of  that 
relationship,  is  derived  as  a  replacement  for  the  proposed  FAE  equation  of  state.  This 
equation  is  not  linear  in  the  forces'  casualty  ratio,  as  is  the  original  FAE  result.  It  describes 
the  ratio  of  the  difference  in  the  square  of  the  forces'  initial  and  final  strengths.  This  new 
relationship  is  consistent  with  both  the  deterministic  and  stochastic  forms  of  Lanchester's 
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Square  Law.  It  is  also  consistent  with  the  historical  data  analysis  due  to  Helmbold  and 
presented  in  Figure  1  above. 

As  a  direct  result  of  the  present  work,  the  July  2005  draft  of  the  TTCP  report  [1]  has  been 
revised  to  include  the  FAE  Equation  of  State  derived  here. 
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